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Abstract 

We calculate the relic abundance of Higgsino-dominant lightest su- 
perparticles, taking account of coannihilations with the superparticles 
which are almost degenerate with the lightest one. We show that their 
relic abundance is reduced drastically by the coannihilation processes 
and hence they are cosmologically of no interest. 



The lightest superparticle (LSP) in supersymmetric (SUSY) standard 
models, mostly taken to be a neutralino, is an intriguing candidate for dark 
matter JI|, 0, ||. An important task is therefore to compute the relic abun- 
dance of the LSPs, which is conveniently represented by Q x h 2 where Q x is 
the ratio of the mass density of the LSPs to the critical one to close the 
Universe and h (0.4 < h < 1) stands for the Hubble constant in units of 
100 km sec^Mpc- 1 . 

This issue has been discussed extensively in the literature ]|, ||, || 0], 
assuming (1) the minimal particle content of the SUSY standard model, (2) 
a GUT relation on gaugino mass parameters (see below), (3) the neutralino 
LSP and (4) R-parity conservation which guarantees the stability of the LSP. 
The abundance of the neutralino LSPs crucially depends not only on their 
mass but also on the composition of them, since the LSP is a linear com- 
bination of neutral gauginos, B, W%, and neutral Higgsinos H\, H 2 , and 
each of them has different pair annihilation processes which determine its 
abundance. 

The results obtained previously are summarized as follows: let us first 
consider the case where the LSP is lighter than the W-boson. If it is almost 
a pure gaugino or a pure Higgsino, its density can reach (or even exceed) 
the critical one and hence it is a good candidate for the dark matter of the 
Universe. If the LSP is an admixture of the gauginos and the Higgsinos, the 
relic abundance is generally small. On the contrary, in the case where the 
LSP is heavier than the W-boson, the annihilation mode to a W-pair opens 
if it is a pure Higgsino or a mixed LSP and the relic abundance is very small 
until the mass increases to a TeV region where Q x h 2 of the LSPs becomes 
again of order unity. If the LSP is a gaugino, there is no annihilation process 



1 



to the W-pair and hence it is still abundant in the Universe. If its mass 
exceeds ~ 500 (GeV), it will overdose the Universe and such a parameter 
region is excluded. 

The previous calculations have, however, included only annihilations 
of the LSP pairs. Griest and Seckel || have pointed out when next-to-the- 
lightest superparticles (NSPs) are slightly heavier than the LSP, this naive 
treatment fails to give a correct answer. Namely, the abundance of the NSPs 
is comparable to that of the LSPs around the freeze-out epoch so that we have 
to take account of annihilation processes involving the NSPs, which they have 
called coannihilations. In ref. ||, they have considered a rather accidental 
case where squarks are degenerate in mass with the LSP, and showed that the 
relic abundance of the LSPs is greatly reduced due to the coannihilations of 
the squarks. We should stress, however, the mass degeneracy rather naturally 
occurs for a Higgsino-dominant LSP |7], |9|. This is because if the electroweak 
gauge symmetry were not broken, all Higgsino states would be degenerate in 
mass. In the broken phase, the mass splitting comes from the mixing with 
gauginos, which are small in the Higgsino-dominant LSP region. Therefore 
the coannihilation effects will be important in the region. 

In this paper, we will evaluate the relic abundance of the neutralino 
LSPs taking account of the coannihilation processes. We consider only the 
processes where neutralinos and charginos pair-annihilate into a fermion pair. 
Since their effects are particularly important for the LSP weighing less than 
the W-boson, we will concentrate on this case. In a region where the Higgsino 
component dominates the LSP, coannihilations of the LSP and the NSPs can 
be mediated by Z- and W-boson exchange. We will show in this region, the 
coannihilations drastically change the previous conclusions: the mass density 
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of them becomes very small, typically Q x h 2 being much less than 10 2 . 
The mass matrix for the four neutralinos is given by 

/ Mi — rriz sin 9w cos (3 mz sin 9w sin (3 \ 

M 2 mz cos 8w cos (3 —mz cos 9w sin (3 

—m z sin 9 W cos /3 m z cos 6^ cos /3 — // 

\ m 2 sin 6^ sin /3 — cos 6w sin /3 — // / 

in the (B, W3, Hi, H 2 ) T basis. Here Mi and M 2 are the mass parameters of 
B and W3, respectively, /i is the supersymmetric Higgs mass parameter and 
tan p represents the ratio of the expectation values of the two Higgs bosons 
[HJ. As is done in most of the literature, we will impose the GUT relation 



for the gaugino mass parameters 

5 

Mi = - tan 2 6 W M 2 ta 0.51M 2 , (2) 
3 

which is derived by solving one-loop renormalization-group equations with 
Mi = M 2 at the GUT scale. As a phase convention, we take M 2 > and 
fi both positive and negative. It is straightforward to diagonalize the mass 
matrix. The lightest of the neutralinos denoted by xi, which is also assumed 
to be the LSP, is 

Xi = Z n B + Z U W 3 + Zi 3 Hi + Z U H 2 , 

E*u 2 = l- (3) 

i=l 

Using Zu, we define the Higgsino purity as 

P = Zil + Z*. (4) 

When p is nearly one, say p >, 0.99, xi is almost a pure Higgsino, which 
is realized if M 2 ^> \p,\. If M 2 the LSP is gaugino-like, whilst for 

M 2 w sa mz it is a general mixture of the four neutralinos. 
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In fig. 1, we have plotted the Higgsino purity of the LSP, p, defined in 
Eq. (H) with tan/3 = 2 fixed. In this paper, we will call the LSP (almost) 
a pure Higgsino if p ^> 0.99, and a Higgsino-dominant LSP if p ^> 0.9. 
The region where 0.99 > P > 0.9 will be called the Higgsino-dominant 
mixed region. Fig. 2 shows the ratio of the mass difference between the LSP 
and the NSP (which is a chargino or a neutralino) to the LSP mass: A = 
(m X2 — m Xl )/m Xl where m Xl and m X2 are the masses of the LSP and NSP, 
respectively. We can see the severe mass- degeneracy in the pure Higgsino 
LSP region. In most of the parameter space, the NSP is the chargino. This 
is also observed when we take a different value for tan f3. 

In calculating the relic density including the effects of the coannihi- 
lations, we use the method developed by Griest and Seckel ||. Here we 
summarize it briefly. Let x% (i — 1> ' ' ' > ^OQ be superparticles with mass m Xi 
(m Xl < ■ ■ ■ < m XN ) and suppose that they are nearly degenerate in mass 
with the LSP, Xi- Due to the mass degeneracy the number density of the 
i-th particle (i > 1), which eventually decays to the LSP, is comparable with 
Hi around the freeze-out epoch. The Boltzmann equation for the total of the 
number densities n = J2i n i is written 

^ = -3Hn - ^(oijVreijinirij - nfnf), (5) 

where is the pair annihilation cross section of the particles x% an d Xji 
v re i is their relative velocity, (• • •) denotes the thermal average and n^ q is 
the number density of Xi i n thermal equilibrium. Since the reactions which 
interchange the superparticles Xi 8 with each other occur much more rapidly 

than their annihilations, the ratio of x% density rij to the total density n is 

1 Note that a particle and its anti-particle should be counted separately when they 
constitute a Dirac fermion. 
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well approximated by its equilibrium value: rii/n « nf 1 /n eq . This greatly 
simplifies the Boltzmann equation as 

^-3Fn-( Wei )(n 2 -Kf). (6) 



Therefore, we can solve eq. (y) by the standard method | | using the effective 
cross section defined by 

<V/ ( 7 ) 

where r, represents the Boltzmann suppression of the density of the heavier 
particle x%- Explicitly 

n eq a-(l + A-) 3/2 e~ AiX 

n e i y ,'1,11 I A ; ) " ^ A ' W 

A i = ( m x« - m xi)/ m xi> ( 9 ) 

where ^ is the degree of freedom of the Xi an d £ = m Xl /T with T being the 
photon temperature. 

The relic abundance of the LSPs at the present day can be calculated 

as 

„ , 2 1.07 x K^GeV- 1 , . 

= — 2 — rTT1 -, (10) 



gj mpij™{a e f f v re i)x 2 dx 



where Xf = m Xl /Tf, Tf is the freeze-out temperature, g* is the effective 
degree of freedom at the freeze-out epoch |Tl| and mpi denotes the Planck 



mass, 1.22 x 10 19 (GeV). In eq. fllpp , we can see the relic density is roughly 
proportional to the inverse of the effective cross section. 

The LSPs freeze out at Xf ~ 20, so when A, (i > 1) is less than about 
0.1 the coannihilation effects are in general important. When the magnitudes 
of the annihilation cross sections involving the heavier particles are similar 
to that of the LSP pair, the coannihilation effects change the relic density 
at most by several factors. On the other hand, if some cross sections with 



the heavier particles are much larger than that of the LSP pair, the relic 
abundance can be reduced by several orders of magnitude. It will turn out 
that this is the case when the LSP is Higgsino-dominant. 

Before giving our numerical result on the relic density, it is illustrative 
to give a crude estimate for the effects of the coannihilations in the degenerate 
limit M2 3> \fi | (see eq. (1)). In this limit, the lightest neutralino is nearly 
either of the following states: 

H SA = ±={H l ±H 2 ). (11) 

The coupling of (neutralino)- (neutralino)- (Z-boson) is given 

h y g/cos9 w )Z ll {H ll ^H 1 - H 2 -fj 5 H 2 ) 

= ^(g/cos6 w )Z fl !l A Yl 5 Hs, (12) 

while that of (neutralino)- (chargino)-(W-boson) is 

9 -W-{H Al »H + + H sl ^H + ) + h.c, (13) 

where g is the SU(2)l gauge coupling constant. Eq. dl2|) implies that in 
a region of the Higgsino-dominant LSP, the coupling of the LSPs to the Z- 
boson is very suppressed: it is proportional to — Z x | ~ m|/(M 2 /i), which 
vanishes at M 2 — > oo limit. 

The annihilation of the LSP pair to fermions is dominated by this sup- 
pressed Z-boson exchange process, since the couplings of the other processes, 
i.e. the Higgs boson exchange and the sfermion exchange, are even smaller. 
Moreover the s-wave annihilation to a fermion pair is suppressed as the initial 
state is a pair of the identical Majorana fermions. Therefore the authors of 
the previous papers have concluded that the Higgsino-dominant LSPs which 
are lighter than the W-boson are (too) rich in the Universe. 
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The reality is, however, they disappear through the coannihilations. To 
see the coannihilation effects, consider the coannihilation process H$,aH ± — > 
(fermions) mediated by the W-boson exchange in the s-channel. The cross 
section for this process is much larger than that of the annihilation of the 
LSP pair discussed above, because there is no suppression factor such as 
Z-J — ZjJ in the couplings to the weak boson and the annihilation occurs in 
the s-wave. Indeed it is estimated as 
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(crv rel ) ~ —"—/((J,) ~ 5/(/i) x 10- b GeV- 2 (14) 

with /(/i) = (/i / m w ) 2 j (4(/x / mw) 2 — l) 2 - The contribution to the effective 
cross section is therefore given by 

(crv rel )r c ~ 5f(jj,)r c x lO^GeV" 2 , (15) 

where rc ~ e~ A ° Xf ~ e _20Ac . Ac and rc represent the mass degeneracy and 
the Boltzmann suppression factor for the chargino, respectively (see eqs. @ 
and (8)). A quite similar estimate can be obtained for the cross section with 
the second lightest neutralino, which is in the same order of the magnitude. 
Since, as we mentioned above, the chargino is the NSP in a large part of the 



parameter space, the effective cross section is dominated by eq. (|I5), which 
yields the estimate of the relic abundance 

2 5 x 10- 10 GeV- 2 IP' 4 f . 

Q x h / v JT^- 16 ) 

{(TeffVrel) r C f(fl) 

Therefore, we can expect that for <, 0.2 (i.e. rc >, 10 -2 ), the relic 
abundance of the LSPs will become very small due to the coannihilations. 

It is tedious but straightforward to calculate the effective cross section 
accurately for a general set of the parameters (/i, M 2 ). In fig. 3, we show the 
abundance of the neutralino LSPs, including the coannihilation processes. 



We have chosen tan (3 = 2, all squark and slepton masses equal to 1 TeV 
and the pseudoscalar mass = 1 TeV. Since the coannihilations occur 
mainly through the Z- or W-exchange, their effects do not depend on the 
choice of the masses of the sfermions and the pseudoscalar. We can obtain 
qualitatively a similar plot for a different value of tan/3. For comparison, 
we have given in fig. 4, the same plot but without taking account of the 
coannihilation processes. By comparing the two figures, the importance of 
the coannihilations is manifest. In the pure Higgsino region with the purity 
p >, 0.99, the mass degeneracy is severe, i.e. A < 0.1, and the coannihi- 
lations greatly reduce the neutralino density. In most of the region, we can 
see fi x h 2 <C 10~ 2 . Thus the LSPs cannot constitute any form of the dark 
matter, i. e. not only the dark matter of the universe but also the dark matter 
of the galactic halos. This contrasts with the conclusion obtained previously 
where the density of the light Higgsinos can reach the critical one. Even 
for the Higgsino dominant-mixed region (0.9 < p <: 0.99) where the mass 
degeneracy is not so severe, the coannihilation effects are significant. They 
reduce the density with an order of magnitude or more. 

In this paper, we have considered the coannihilation effects to calculate 
the present abundance of the neutralino LSPs whose mass is less than that 
of the W-boson. We have shown that in the Higgsino-dominant LSP region 
the relic density gets very small. On the other hand, as mentioned above for 
a Higgsino-dominant LSP heavier than the W-boson, the annihilation to the 
W-pair is kinematically allowed and fi x h 2 is much smaller than order unity 
until the LSP mass reaches the TeV region. Combining these two results, 
we can conclude that the Higgsino-dominant LSP is no longer an interesting 
candidate for the dark matter. 
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We would like to make two comments: firstly our calculation does not 
contain annihilation processes whose final state is other than a fermion pair. 
When Higgs bosons are light, the LSPs and NSPs will annihilate to them sig- 
nificantly. If this is the case, the relic density becomes even smaller. Moreover 
we may consider a process such as 

H SA H ± ^W ± ^^W ± , (17) 

which would enhance the effective cross section by several factors. Although 
it will not change our result drastically, a more detailed and accurate analysis 
is significant |fl2| . Secondly there are other cases where the coannihilations 



are important. For example, when the GUT relation for the gaugino masses 
(|2J) is relaxed, the neutral Wino can become the LSP Hj. Since it is highly 
degenerate in mass with its charged counterparts, the coannihilations involv- 
ing them dominate the effective annihilation cross section and hence their 
relic abundance is very small. This issue will be discussed in a separate 
publication [jlTfl . 

We would like to thank T. Yanagida for careful reading of the manuscript 
and useful comments. We are also grateful to H. Murayama and T. Yanagida 
for discussions. 
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Figure captions 

Fig. 1. The Higgsino purity of the LSP, p , defined by p = Zf 3 +Zf A ; (a)/i > 
and (b)/i < 0. We have chosen tan/3 = 2. "LEP" means the excluded 
region by LEP constraints and "M x > M^" means the region where 
the lightest neutralino is heavier than W-boson with which we are not 
concerned. 

Fig. 2. The ratio of the mass difference between the LSP and the NSP to 
the LSP mass, A , defined by A = (m X2 — m Xl )/m Xl ; (a) fi > and 
(b) fi < 0. We have chosen tan/3 = 2. The meanings of "LEP" and 
"M x > M\y" are the same as in fig. [I]. 

Fig. 3. The relic abundance of the neutralino LSPs, including the coanni- 
hilation processes; (a) \i > and (b) fi < 0. We have chosen tan/3 = 2, 
all squark and slepton masses equal to ITeV and the pseudoscalar mass 
itla = ITeV. Note that the coannihilation effects are independent of the 
squark and slepton masses and the pseudoscalar mass. The meanings 
of "LEP" and "M x > M^" are the same as in fig. |I[ By comparing 
this with fig. 4, we can find that the coannihilations greatly reduce 
the relic abundance of the LSPs in the Higgsino-dominant LSP region 
(M 2 > |/i|). 

Fig. 4. Same as fig. |3] but without taking account of the coannihilation 

processes. 
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